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Abstract 

We study the sensitivity of the Tevatron and the 7 TeV LHC to a leptoquark S coupling to a top quark and a 
charged lepton L (= e, fi, or r). For the Tevatron, we focus on the case ms < nit, where the leptoquark pair 
production cross section is large, and the decay is three-body: S — > WbL ± . We argue that existing Tevatron 
observations could exclude ms 160 GeV. For mg > m t , we show that the LHC experiments with low integrated 

> , luminosity could be sensitive to such leptoquarks decaying to tl with I = u or r. 
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, 1 Introduction and Review 

Leptoquarks [1] are bosons which couple to a lepton and a quark. Although they are not known to address current 
issues in particle physics (such as the identity of dark matter or the hierarchy problem), they can be motivated in 
several ways. Most pragmatically, they are strongly interacting and their decay products include leptons, so they are 
interesting search candidates for hadron colliders. The Tevatron sets bounds on leptoquarks which decay to first and 
second generation fermions, and to 6s; in this note, we consider leptoquarks which couple to the top quark and any 
£X| charged lepton L ± (L s {e, n, r}). We discuss the bounds that could be set with 4.3 of Tevatron data, and the 

> ! prospects for the 7 TeV LHC with 1 fb~ x . 
Leptoquarks can arise in several extensions of the Standard Model, such as Grand Unified Theories [2], Technicolour 
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[3] and i?-parity violating Supersymmetry[4]. We focus on scalar leptoquarks called S, with baryon and lepton number 
conserving interactions, and a mass m$ 1 TeV. Several recent models [5, 6, 7] include such leptoquarks. The 
Lagrangian describing their renormalisable interactions with Standard Model (SM) fermions and singlet neutrinos v 
is [8] 

£lq = S {\ Ll s o iiT2q c + A RSo eu c ) + S X R ^ o ed c + {\ls 2 1u + ARs 2 eg[zT 2 ])S , 2 + A L § 2 M? 2 + l[iT 2 ]fq c ■ S 3 
> +S' \^ So Vu c + Sft'zVd + \^ S2 Vq\iT 2 }S 2 + h.c. (1) 

• rH ' 

■ where the leptoquark subscript is its SU(2) representation, the As are 3x3 matrices in the lepton and quark flavour 
^ | spaces and are labelled by the SU(2) representation of the leptons [L = doublet, R — singlet) and the leptoquark 
name, r 2 is a Pauli matrix (so ir 2 provides the antisymmetric SU(2) contraction), the SM SU(2) singlets are e, u, d and 
and in this equation, q and £ are the doublets. For most of the rest of the paper, L and £ label physical particles 

£e{e,fi}, i£{e,/i,r} 

In eqn (1), we included for completeness, leptoquarks which couple to singlet neutrinos vr. If the neutrino masses are 
Dirac, these interactions could allow S — > tv without S — > bv. However, we do not analyse such decays. Notice that 
we neglect, or set to zero, the (renormalisable) interactions of the leptoquark with the Higgs, which naturally should 
be present, and can contribute via loops to precision electroweak parameters [9] and neutrino masses [10]. 

To look for leptoquarks, some theoretical expectations about the structure and hierarchy of their interactions would 
be helpful. Various theoretical arguments can suggest that the largest leptoquark couplings should be to the third 
generation, at least in the quark sector. This arises, for instance, in the Cheng-Sher ansatz [11] for flavoured couplings 



»"«,/^ (2) 



where v = 175 GeV is the Higgs vacuum expectation value. This would give the largest leptoquark coupling to t and 
t, and can arise Randall-Sundrum type extra dimensional models, or in composite models as recently discussed in [6]. 
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Figure 1: Lowest order diagrams for leptoquark single and pair production; single production can be neglected for a 
leptoquark which couples only to t quarks, because the t content of the proton is small. 

The expectations of this ansatz are compared to current low energy constraints in [12]; improving the sensitivity of 
K — > nvv could probe this pattern for leptoquarks that couple to neutrinos. 

A phcnomenological "bottom-up" approach to the couplings of new particles, motivated by the success of the 
Minimal Flavour Violation hypothesis [13, 14], is to construct them by multiplying the known mass and mixing 
matrices of leptons and quarks. Some possibilities for leptoquarks were studied in [15]. In this approach, Nikolidakis 
and Smith [16] noted that a New Physics coupling with a single lepton index L can be proportional to 

e LJK [Y e Y}m v ] JK (3) 

where Y e is the charged lepton Yukawa matrix (index order doublet-singlet), and m v is the majorana neutrino mass 
matrix. This idea was studied for leptoquarks in [15]. Since m v is fairly democratic, the Y e Y^ hierarchy selects the 
t index. The totally antisymmetric SU(2) tensor is e, so this construction favours couplings to the e and /i. It is 
therefore interesting to study leptoquarks which decay to t and any charged lepton: L = r, jj,, or e. The /i and e are 
particularily attractive final state particles for hadron colliders: if one of the Ws from the ts decays to the I = e or //, 
as occurs <~ 29 % of the time, the final state would be jets +|?t + ^ ± ^ ± ^ T (see figure 2.) 





Figure 2: Possible decay chain for a pair of scalar leptoquarks interacting with tops and muons. 

There are various experimental constraints on leptoquarks. At hadron colliders, they can be singly or pair produced 
via their strong interactions (see figure 1). As discussed in [17], single production can lead to the same final state as 
pair production. However, since we are interested in leptoquarks that couple to the top quark, we can neglect single 
production, because the top density in the proton is negligeable. The cross section, for pair production from gg or qq, 
has been computed at Next to Leading Order (NLO) [18], and included in the PROSPINO program, which we used to 
produce figure 3. The Tevatron [19] has searched for leptoquarks decaying to any lepton and a quark other than the 



top, with a coupling A <; 1CP 8 . The restriction on A ensures that the leptoquark decays at the collision point. The 
bounds depend on the final state; a recent review [20] gives ms *S 210[r6], 214[j/<7], 247[^6], 299[eg, eb], 316[^tg, /ib] 
GeV, where q € {u, d, s, c}. Leptoquarks have also been searched for at the HERA ep collider, which allow to exclude 
a s-channel resonance with A 0.1 and ms 250 — 300 GeV [21]. Finally, there are bounds on two quark, two lepton 
contact interactions from several (mostly accelerator) experiments [22], which give interesting constraints (see e.g. 
[12], [23]) on leptoquarks interacting with lower generation fermions. The prospects of discovering leptoquarks above 
the various backgrounds at the early LHC have been discussed in [24]. 

There are numerous precision/rare decay bounds on leptoquarks, which usually apply to products of different As, 
and depend on the SU(2) representation of the leptoquark. Some recent compilations are [25] (bounds from meson 
anti-meson mixing, allowing for complex couplings) and [12] (mostly tree processes). In general, it is clear that these 
bounds exclude flavour-democratic \ L ® ~ 0(1) for ms < TeV. Certain processes, such as K — > ■nvv or Kl — > /i ± e =F 
provide much more stringent bounds. 

Bounds and prospects for a "third generation" leptoquark interacting with a top, have been discussed by several 
people, in particular Eboli and collaborators. The constraints from the loop contribution to leptonic Z decay [26] are 
satisfied if A e for ms ~ 300 GeV (for both the leptoquarks of eqn (1) which couple to tn). To extrapolate this 
bound for leptoquarks in the range 300 GeV — > my/, we assume that the bound can be scaled as X/ms ^ e/ (300GeV), 
see eqn (4). The LHC prospects of a leptoquark decaying to tr or br were discussed in [27], who emphasized the 
interesting one, two and three lepton final states which could be detected above backgrounds. Gripaios et al recently 
implemented the various leptoquarks of eqn(l) in herwig [28], and discussed kinematic reconstruction techniques for 
leptoquarks decaying to third generation fermions at the 7 TeV LHC. 

In this paper, we study leptoquarks which couple to tops, but not to 6s or lower generation quarks, because 
leptoquarks with an 0(1) branching ratio to bv, or jet + e or /i, are already excluded by the Tevatron up to ms ^ 
200 — 300 GeV [19, 20]. We are therefore interested in leptoquarks which couple to singlet up- type quarks, that is, the 
SU(2) singlet leptoquark So with coupling Xrs (and \ls = 0), or the doublet leptoquark S2 with couplings Xls 2 7^ 
and Xrs 2 = 0. Neither of these leptoquarks arise in i?-parity violating Supersymmetry. We then restrict the coupling 
to third generation quarks, and assume, for the body of the paper, a branching ratio of 1 to the final state of top + 
the charged lepton L^ 1 of our choice. 

The NLO cross section for leptoquark pair production [18], via the strong interaction, is plotted in figure 3. This 
shows that the Tevatron with 5 fb^ 1 of data could produce <; 0(5000) pairs of leptoquarks with ms Ss 150 GeV, 
whereas the 7 TeV LHC with 1/6 -1 could produce a thousand pairs of 300 GeV leptoquarks. Section 2 outlines a 
simple counting experiment that compares current Tevatron data, to leptoquarks with ms < m t , decaying via an 
off-shell top to bW and a charged lepton. It suggests that the Tevatron could exclude such leptoquarks, for leptoquark 
masses sufficiently below m t . In section 3, we briefly mention prospectives with 1 of data from the 7 TeV LHC. 
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Figure 3: Pair production cross section for SU(2) singlet leptoquarks at the Tevatron and the 7 TeV LHC. Standard 
Model ti production, with <J pp ^tt(Vs = 7 TeV) ~ 165 pb and (J p p^tt{\fs = 1-96 TeV) ~ 7.8 pb, could be a significant 
background to the leptoquark signal, in particular in the mg ~ to* ± 15 GeV region, where the leptons leptoquark 
decays would be soft. 



2 Leptoquarks with ms < m t at the Tevatron 



In this section, we consider leptoquarks S interacting only with the top (and any L), with masses in the range 
m\y + ni}, < ms < 2mw < m t . They would be copiously pair-produced via their strong interactions at the Tevatron. 
There are two reasons for this limited mass range, despite that figure (3) suggests several hundred lcptoquark pairs 
could be produced at the Tevatron up to masses ms ~ 250 GeV. Firstly, in the range mj ~ m t ± 15 GeV, the lepton 
produced with the almost-on-shell top is unlikely to pass the pr > 15 — 20 GeV cuts that we impose. Secondly, it is 
convenient to analyse separately the ms < m t and ms > m t cases; so we study the former at the Tevatron and the 
latter at the LHC. 

Leptoquarks with ms < m t could also be singly produced in top decay; however we neglect this process, because 
BR{t — > SL) < 2|A| 2 /|y t | 2 is suppressed by the leptoquark coupling A (y t is the top quark Yukawa coupling). We 
consider the range 

where the upper bound is approximately the constraint from leptoquark loop contributions to the ZLL vertex[26]. 
It implies that BR(t — »■ SL) is negligeable, when it is kinematically allowed. The lower bound ensures that A is 
sufficiently large that the leptoquarks decay at the collision point (no displaced vertex); see the discussion after eqn 
(14). We also assume that S interacts only with a top and a charged lepton L ± , so BR(S — > bWL) = 1, where L is a 
e, or t. 

Since the leptoquark decays to a singlet tn, its decay rate to bWL, via an off-shell t, has a simple analytic form, 
which we obtain in subsection 2.1. This allows us to implement the three-body decay in PYTHIA, as the product of 
the two body decays T(S — > t*L)T(t* — > bW) with variable m t and leptoquark coupling A. This is discussed at the 
end of section 2.1. 

Subsection 2.2 contains preliminary estimates of the contribution of such leptoquarks to the jets + Ifix + t [29] 
and jets + IpT + ^t^J [30] data sets used by DO to measure the ti production cross section (I here means e ot fx). We 
consider separately the cases S — > tr ± and S — > i^ ± ; we assume that the bounds which could be obtained on5-^ te ± 
are similar to those on S — > tjj^. 

2.1 The decay rate S — > bWL^ 1 for m s < m t 

If the masses of the b and L E {e, /i, r} are neglected, then the invariant mass of the bW system (or equivalcntly, the 
magnitude of the four-momentum carried by the off-shell top in the decay S — > bWL), is 

f2 = m lw = (Pb + Pw) 2 = %Pb ■ Pw + m 2 w . (5) 
The differential three-body decay rate can be written [31] 



~dt 



= (^16^|/^| 2 ^l^l^ ^ 



where the L parameters are in the S rest frame, and the starred b parameters in the bW rest frame. 
The matrix element for S — > bWL is 



(where ul is the spinor field for L), and is simple in squared form because the top must flip chirality on the internal 
line: 



l .,,2 m 2 t \ 2 g 2 ( p L -pwPb-pw \ fa . 



To evaluate the angular integrations of eqn (6) with |.M| 2 from eqn (8), requires the Lorentz transformation of the 
b 4-momentum in the S frame (pb), to the t frame(pj). Writing 

E b = 7 £ 6 *(l+/3cos0*) (9) 



A 



with 7 = E t /m t , (3 = \p t \/m t , and using 



1 f-> m w\ i- i m S 



gives 



I^T I — 2 ( 1 ^ ) , \PL\ = ^[1- — ) , (10) 



n l». + 2 PL - P Z?- PW ) d%<m L = 4*>m% (l - ^ ) — ( 1 - 1^ ) ( I + 2^- , 



2 

"*W / V m S/ ""W 

So one obtains 

dT T(S -> t*L) r(f -> W6) 



df 2m t 7rm t (t 2 - m 2 ) 2 + m 2 T 2 y ' 

where the two body decay rates of the lcptoquark S and the top quark t are the rest-frame formulae, with rrit replaced 
by t: 

nt^Wb) = ^^(l- n f-) 2 (l + 2^-) , T(S^ t *L) = ^(l-*Y . (13) 



64tt m 2 w \ t 2 J \ t 2 J v ' 16tt 

A check that can be performed on eqn(12) is to take the limit t 2 — > mf. Using the identity: 

7T x + (r 

with x = t 2 —m 2 , the integration can be performed, and one obtains the two-body leptoquark decay rate T(S — > tL), 
as expected. We are interested in the case t 2 — m 2 3> T t mt, so we drop the T t term in the denominator of eqn (12). 
The total decay rate T(S — > bWL) can also be obtained analytically, but is not illuminating. It is plotted on the left 
in figure 4 for A = 1. The leptoquark will decay in less than a centimetre for A <; 10~ 3 at ms — 100 GeV, and for 
A > few x 1CT 6 at m s ~ 160 GeV. 




Figure 4: On the left, the total decay rate T(S — > bWL), for A = 1, as a function of the leptoquark mass ms, where L 
is a charged lepton e,/x or r. On the right, the fraction of decays to a bW of invariant mass t for various leptoquark 
masses ms- 

Equation (12) implies that the decay S bWL can be computed as the two-body decay S —> t*L, where t* is a 
top quark of mass t, followed by the decay t* — > bW, provided the whole process has a i-dependent coupling constant 
oc l/(t 2 — m 2 ) 2 . We implement the ms < m t leptoquarks in pythia from this perspective: for each event, the top 
mass is randomly selected, distributed between mw + mb and ms according to eqn (12). The two-body leptoquark 
decay is then performed by pythia followed by a two-body top decay. We checked that the resulting bW invariant 
mass distribution reproduces eqn(12). 



2.2 Potential bounds from the Tevatron 



In this section we estimate the number of events from m s < m t leptoquarks, that could appear in two recent DO data 
samples used to measure the ti production cross section: ^ ± + ]fij< + jets [29], for leptoquarks decaying to tr ± , and 
tftj + $t + jets [30] for leptoquarks decaying to . These two analyses arc not optimal to select these leptoquark 
signals because they require exactly one and two leptons, respectively, in the final state. However, due to the important 
leptoquark pair production cross section when ms < m t , we will show that a significant number of leptoquark events 
would nevertheless enter in these data samples. 

We obtain the leptoquark events using PYTHIA6.4 [32] with tauola [33] to decay rs 1 , and reconstruct jets using 
the anti-fci algorithm of the fastjet package [34] . All final state particles except neutrinos and charged leptons are 
used to construct the jets. We do not include a detector simulation. We only impose the preselection cuts of the 
experimental analyses on our leptoquark events, and count the number of remaining events. DO uses multivariate 
techniques to further discriminate ti signal from W+ jets or Z+ jets, which we do not consider. We study two cases: 
S -» tr~, and S -» tfi + . 

In pythia, leptoquarks can decay to t and L + , but not 2 to t and L~ . The first is appropriate to the SU(2) doublet 
leptoquark S2, whereas the second would correspond to the singlet leptoquark So- In the S — > t\i case, the events 
corresponding to tyT or t/i+ should be equally detectable, so we let pythia decay the leptoquarks to tfi + , and apply 
the resulting bounds to both S2 and So- 

The case of S — > tr^ is more delicate, because the angular and energy distribution of the tau decay products 
depends on the charge and polarisation of the r. We therefore ask TAUOLA to flip the sign of rs produced in So decays. 
In the case of leptonic r decays, energetic l~ are emitted preferentially anti-parrallel to the direction of motion of 
a relativistic tr (see e.g. the e distribution from ji decay in [31]). So for the singlet leptoquark So, which decays 
to tji(Tji)~ , the neutrinos from the leptonic tau decays frequently carry most of the px of the r. Since the rs are 
already not very energetic, this means the charged leptons from their decay may not pass px cuts, so hadronic r decays 
are more useful. We assume that the bounds we obtain on the singlet leptoquark So can be applied to the doublet 
leptoquark S2, which decays S2 — >■ t(rx) + , because our bounds will mostly come from events with hadronic r decays 3 . 

2.2.1 S-J-tr* 

Consider first the production of a leptoquark anti-leptoquark pair, followed by leptoquark decay to a t and a t~ or t + . 
This could contribute to the £ ± + Ipx + jets signal from which DO extracted the Standard Model top pair production 
cross section cr t j in [29]; since the leptoquark process should have more jets than ti production, we focus on the £ ± and 
at least 4 jet sample. 

In our simulated sample of leptoquark events, we require that a W ± from the t or t decays to v and e ± or fj,^, 
which reduces the cross section by a factor 

BR(£ ± +$ T + n jets) = .22{W+ -> t+v) x M(W~ -> had) x 2(W+ W~) a .29 (15) 

Then we impose the following cuts, patterned on the preselection of the DO analysis. We require 

1. IpT > 25 GeV 

2. a lepton with px > 20 GeV and \t]\ < 1.1(e), 2.0(h) an d no second lepton with px > 15 GeV 

3. at least 4 jets with p T > 20 GeV, and \rj\ < 2.5 

With e S i m the fraction of simulated events which pass these cuts, the inclusive leptoquark signal efficiency is simply 
e(ftr, It, 4j) — Esim x BR. This efficiency is given in colomn 3 of Table 1. We then estimate the number of leptoquark 
induced events in the 4.3 of data used in [29] to be the last colomn of table 1. 

The total number of observed [expected] £ ± + $t+ > 4 jets events in [29] is 1795[1796± 158]. Using the modified 
frequentist CLs method [35] , the number of signal events excluded at 95% C.L. is 388. By interpolating our results, 
we find that leptoquarks with mass 

m s < 158 GeV for BR(S tr ± ) = 1 

are excluded at 95% C.L. 

1 We modified the TAUOLA-PYTHIA interface so that it finds and assigns polarisation to rs from leptoquark decay. 

2 The various scalar leptoquarks of eqn (1) have recently been included in HERWIG [28], which would avoid this limitation that PYTHIA 
only knows one chiral structure for leptoquark interactions. 

3 We checked that changing the r polarisation makes a relatively insignificant change to the bounds. 



m s (GeV) 


a (pb) 


e(#r,lMj) 


N(LQ) 


160 


1 


.0823 


367 


140 


2.4 


.0618 


658 


120 


6 


.0389 


1035 


100 


16 


.0149 


1060 



Table 1: The second colomn is the leptoquark pair production cross section at the Tevatron for various masses. The 
leptoquarks decay to tr^. The third colomn estimates the fraction of events remaining after the cuts given in section 
2.2.1. The last colomn is the expected number of leptoquark-induced events in the DO lepton + > 4 jets sample [29] 
based on 4.3 



2.2.2 S-^tfj^ 

Consider now a pair of leptoquarks which decay to a t and a ^ , which could contribute to the ifij + $t + jets data 
sample from which DO extracts a t i [30]. In this analysis [30], can be e + e~,e ± fT F or n + p~ ■ The leptoquarks we 

study would contribute mostly to e ± /i T or n + ■ However, to be conservative, we compare our expectations to the 
observed number of if £ J + jets events, including e + e~. 

We simulate leptoquark pair production, followed by leptoquark decays S — > W + b[i + , and require that one W 
decay to a charged lepton e, p, or r, which should represent a fraction 

.34(VK+ -> i + v) x .66(W~ had) x 2(W + W~) ~ .45 (16) 

of the events. We require a leptonic W to ensure missing transverse energy in the event, but include also the W ' — » tv 
decays, because a lepton from the W is not neccessary since two charged leptons are already coming directly from the 
leptoquark decay. 

Our cuts to select two charged leptons are patterned on the DO analysis [30]. We therefore require exactly two 
opposite sign leptons of pt > 15 GeV with \t)\ < 1.1(e), 2. 0(/z). Then, D0 uses a multivariate technique (Bayesian 
decision tree or BDT) to further discriminate top pair events from i?/7*+jets events which we can not take into 
account. But since the topology of our leptoquark signal is very close from the one of top pair production, we believe 
that leptoquark events will nevertheless pass the selection cut applied on this BDT output with an efficiency very 
close to the one from top pair events. In the following, we will not take into account the efficiency of the BDT, and 
simply replace it by a cut Jf,T > 25 GeV. Then, we count the number of jets satisfaying pt > 20 GeV and r\ < 2.5 and 
require at least 3 jets. 

In the third and fourth colomns of table 2, we give the estimated fraction of leptoquark events which would pass 
the above cuts (obtained by multiplying eqn (16), and the fraction of simulated events which pass cuts) and the 
expected number of leptoquark-induced events in 4.3 /6 _1 . These numbers can be compared to the ~ 51[65 ± 15] 
observed [expected] events 4 , bearing in mind that we have not simulated detector effets and that we did not take into 
account the efficiency of the selection cut on the BDT output rejecting Drell-Yan events, which is around 70%. From 
those events, we computed that the number of signal events excluded at 95% C.L. is 39. We therefore see that the 
expected number of leptoquarks events is much larger than this number: this leptoquark signal would significantly 
contribute to the number of events observed in the D0 analysis, and we can conclude that leptoquark with mass 

m s < 160 GeV for BR(S ->• t^) = 1 

are excluded at 95% C.L. 

3 At the 7 TeV LHC 

Leptoquarks decaying to a top and an e or fi are attractive search candidates for the early LHC because the final state 
contains leptons and many jets. If a W from the t or i decays leptonically, various combinations of same sign and 
opposite sign leptons of different flavour can be obtained (see figure 2). 

Since the events contain many jets, the leptoquark pair production and decay should be calculated at NLO, so that 
the Monte Carlo simulation matches as well as possible to the real events. In addition, detailed study of backgrounds 

4 The numbers are extracted from a histogram. 
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m s (GeV) 


a (pb) 


e(# T ,20S£, 3j) 


N(LQ) 


160 


1 


.0900 


387 


140 


2.4 


.0752 


776 


120 


6 


.0500 


1288 


100 


16 


.0090 


960 



Table 2: The second colomn is the leptoquark pair production cross section for various masses. The leptoquarks decay 
to t/i . The third colomn estimates the fraction of events remaining after the cuts given in section 2.2.1. The last 
colomn is the expected number of leptoquark-induced events in the DO £ + £~ + > 3 jets sample [30] of 4.3 



would be required to identify suitable cuts to select leptoquark events and identify the leptoquark mass. We leave 
this analysis to the experimental collaborations, and here, we merely estimate the fraction of events at the LHC, with 
cuts similar to recent LHC tt results [36, 37]. 

We consider leptoquarks that decay with a branching ratio of one to either tr~ , or tfi + , with a mass in the range 
200 — 400 GeV (so they decay to an on-shell top). The production and decay are calculated by PYTHIA 6.4[32]. Jets 
are reconstructed with the anti-fct algorithm of the fastjet package[34], with R = .5. To estimate a total number of 
surviving events, we assume 1 of data. 

3.1 Counting events: S — > tfi + 

Consider first the decay S — > tp + . If one W decays leptonically, this could be searched for in events with £ + £~ + $t+ 
jets. CMS recently determined the tt production cross section [36] (with two leptonic Ws) from such events, and our 
cuts are patterned on this analysis. We expect more leptons and jets in SS production than in tt production, so we 
impose: 

1. $ T > 30 GeV 

2. Exactly two opposite sign charged leptons (e ± ,^i ± ), with pr > 20, \r]\ < 2.5. 
Or alternatively, two OS leptons, with at least one other lepton. 

3. at least four jets of pr > 30 GeV, and \r]\ < 2.5 

The CMS analysis has an isolation cut for the leptons; we instead require that the simulated leptons who pass pr 
cuts be produced in W or S decays (to avoid high px leptons from meson decays). We allow all decays to our Ws 
in pythia. This means, for instance that our simulation now includes events with two leptonic Ws, which could 
pass cuts if there are additional QCD jets (this accounts for ~ 10% of our events at ms = 200 GeV). The fraction 
of events that survive cuts 1, 3 and either of the versions of 2, are respectively defined as e(]pT , = 20S£,4j) and 
e(lpT, > 20S£, 4j), and are given in colomns three and five of table 3. The number of events at the 7 TeV LHC with 
C = l/fr -1 of integrated luminosity is estimated in the fourth and sixth colomns. 



m s 


Cprod/pb 


e($ T ,= 20S£,4j) 


N = (LQ) 


e{jpr,>20S£,4j) 


N>(LQ) 


200 


12.5 


.055 


683 


.035 


438 


250 


3.69 


.095 


352 


.094 


346 


300 


1.3 


.104 


136 


.116 


151 


350 


0.515 


.109 


56 


.12 


62 


400 


0.224 


.121 


27 


.129 


29 



Table 3: The second colomn is the leptquark pair production cross section at the 7 TeV LHC. The leptoquarks decay 
to i/i*. The third colomn is the fraction of events which pass the cuts of section 3.1 with exactly a pair of opposite 
sign leptons, and the fourth colomn is the estimated number of events in 1 fb^ 1 of data. Colomns five and six are the 
same, for the 3 or more lepton cut of section 3.1. 

We can compare to the CMS determination [36] of the tt production cross section, based on 3.1 pb~ x of data from 
the 7 TeV LHC. In the 20S£ and > 4 jet bin, CMS observes one event, where ~ .75 — > 1.5 signal events are expected. 



« 



From [38], it appears that the background is ^ 1/3 of the signal. We anticipate that a 200 GeV leptoquark would 
contribute <~ 2 events in the > 4j and exactly 2 OS lepton bin. 

The integrated luminosity available now (winter 2011) is of order 35 pb~ x , or ten times that used in the CMS analysis 
[36]. This suggests that a t i measurements at the LHC are already sensitive to leptoquarks S with BR(S — > tfi^) = 1 
and ms ^ 200 GeV. Furthermore, searching for > 3£ and > 4 jets would be more sensitive to such leptoquarks. 

3.2 Counting events: S — > tr~ 

Consider now the decay S — > tr~ with BR(S — > tr~) = 1. This decay would be more challenging to reconstruct, 
because the neutrinos from both rs and a leptonic W can contribute to JpT ■ We decay S — > tr + in pythia, and 
tell TAUOLA to flip the sign of the ts from the leptoquarks : r ± — > r T , with helicity assigned as if it were chiral 
singlet (tr). Similarly to the discussion of 5* — > tr at the Tevatron (see section 2.2.1), we attempt to constrain these 
leptoquarks at the LHC from lepton + jets +Ipr events. However, in our simulation of S — > tr at the LHC, unlike 
that of section 2.2.1, we allow all decays to the Ws from the tops. This is because, at the LHC, leptons produced in 
r decay can be energetic enough to pass pt cuts. 

We then impose the following cuts, patterned on an ATLAS [37] analysis which extracts the ti production cross 
section from lepton + jets +$t events: 

1. $t > 25 GeV, where all the neutrinos are summed into fx 

2. at least four jets with > 25 GeV, |?7| < 2.5. 

3. we require at least one e ± ,p, ± with p T > 20 GeV, \rj\ < 2.5. To mimic an isolation cut on our simulated leptons, 
we then check that the e + ,e~,p + and pT with highest ptS are separated from the jets which pass cuts, by 
AR = yj (rjj — r\i) 2 + (<pj — <\>i) 2 > .3. The leptons failing this check are rejected. Then, we may in addition 
require either exactly one lepton (as in the ATLAS analysis), or, at least 7 jets and/or and/or ^ which pass 
these cuts. 

To take into account both hadronic and leptonic t decays, the requirement of > 7 jets and/or leptons is applied. 

Our estimates of the fraction of leptoquark events that would pass these cuts (for the three possible lepton cuts), 
and the number of events in 1 fb~ x of data, are in table 4. Notice that the efficiencies, for finding an S — )■ tr^ 1 
leptoquark in such single lepton events, are higher than the efficiencies to find S — > t/i ± leptoquarks in the dilepton 
events, given in table 3. This is because there is always $t (^vs) in the S — > tr ± final state, and there is an or p ± 
approximately 2/3 of the time. Whereas requiring JpT in the S — > tp ± decays, imposes a leptonic W, which occurs 
~ 29% of the time. 



m s 


cr/pb 


e(> U) 


N(> U) 


e(l £) 


N(l I) 


e(> 71 + j) 


N(> 7£ + j) 


200 


12.5 


.160 


2000 


.143 


1788 


.039 


488 


250 


3.69 


.297 


1096 


.241 


889 


.126 


465 


300 


1.30 


.374 


486 


.285 


370 


.199 


259 


350 


.515 


.428 


220 


.322 


166 


.234 


121 


400 


.224 


.451 


101 


.328 


74 


.264 


59 



Table 4: The second colomn is the leptoquark pair production cross section at the 7 TeV LHC. The leptoquark decays 
to tT^. The third, fifth and seventh are the fraction of events which survive the cuts of section 3.2; the three es 
correspond to the three different lepton cuts. Then the fourth, sixth and eighth colomns are the estimated number of 
events passing cuts in 1 fb~ x of data, for BR(S — > tr~) = 1. 

It is less straightforward to anticipate LHC sensitivity in this channel. ATLAS [37] measured the ti production 
cross section in the single lepton + jets +$t channel, with 2.8 pb~ x of data. In the bin containing > 4 jets of which 
two are tagged as 6s, they observe 37 events, expected 30 from ti, and estimate the background as 12.2 ± 3.9. The 
ATLAS b tagging efficiency varies (40 — 60% for_Ej : 25 — > 85GeV); if we assume that ~ 50% of the 6s from leptoquarks 
are tagged, one can guess from table 4 that a ms — 200 GeV leptoquark, with BR(S — > tr~) = 1, could contribute 
^1 — 2 events to this bin. Explicitely counting events with extra jets (beyond the four expected from ti), or events 
with extra leptons (ATLAS required one and only one), could improve the sensitivity to leptoquarks decaying to tr~. 
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4 Summary 



Like the Higgs, a scalar leptoquark is a boson that couples to two fermions. Since the quark Yukawa couplings to 
the Higgs are hierarchical, and flavour physics in the quark sector follows Standard Model expectations, one can 
anticipate that leptoquarks, like the Higgs, interact preferentially with third generation quarks. This also arises in 
several models. However, expectations for leptoquark couplings to leptons are less straightforward to extract from 
lcpton mass matrices. The charged leptons are hierarchical, so one could imagine that leptoquarks should preferentially 
decay to ir^. On the other hand, the neutrino sector is comparatively democratic, suggesting that leptoquarks could 
decay to t and any lepton. The te^ and tfi^ final states could be interesting search channels for the early LHC. 

This paper studied possible bounds on leptoquarks S, with a mass in the range 100 GeV < ms < 400 GeV, which 
are pair-produced via their strong interactions, and decay to a top quark (and only the top; no b, c...), and a charged 
lepton L = e,/i or r. We expect the Tevatron, with its high luminosity, to be sensitive to the range ms < mt, where 
the leptoquarks decay to the three-body final state bW ± L T . Leptoquarks with ms > m t could be found at the LHC. 
The range ms — m t ± 15 GeV appears difficult: the soft leptons in the final state may not pass px cuts, so the SS 
final state becomes difficult to distinguish from tt. Recall that at ms ~ m t , the SS production cross section is <~ 1/10 
of the tt production cross section. 

We estimated the number of leptoquark-induced events containing lepton(s) plus jets, using pythia6.4, tauola, 
and the anti-fc t jet algorithm of fastjet. We include no detector simulation. We consider separately the sensitivities 
to leptoquarks which decay to tr ± , or i/i ± , assuming a branching ratio of 1 in each case. We further assume that our 
estimates for the t^ final state could apply to leptoquarks decaying to te ± . 

Our results suggest that current determinations of the tt production cross section, both from the Tevatron and the 
LHC, could constrain a leptoquark with BR(S — > i/i ± ) = 1, and a mass of order 100 — > 250 GeV. At the Tevatron, 
DO determines the tt production cross section a t i from the final state l i £j + $t+ > 3 jets, using 4.3 fb^ 1 of data. 
We estimate that these results could exclude 

m s ^ 160 GeV for BR(S^t£ ± ) = l, £ e {efi} 

At the LHC, CMS obtained a t f from events containing iftf + $t+ > 2 jets. It is possible that leptoquarks with 
ms ^ 200 GeV, could have contributed a few events to the scantily populated (1 event) > 4 jet bin. 

Leptoquarks with BR(S — > tr ) would contribute to the the final state i + Jpx + > 4 jets (which is used to 
determine a t i), if at least one of the Ws or rs decays leptonically. A significant fraction of the leptoquark events 
should have more than four jets, but the available data sets present a single > 4 jet bin, so we cannot profit from this 
property. We estimate that a DO analysis could exclude 

m s £ 160 GeV for BR(S -> tr*) = 1 

ATLAS also obtained u t i from events with £ + $t+ > 4 jets, but leptoquarks with mj ~ 200 GeV would be consistent 
with the backgrounds. 

The current integrated luminosity of the LHC is significantly larger than that used in the analyses we compared 
to [37, 36]. So we anticipate that the winter 2011 determinations of u t i at the LHC should have some sensitivity to 
the leptoquarks discussed here. However, SS production followed by S — > tL^, should usually give a final state with 
more leptons and/or jets than tt production. This means that analyses counting additional leptons and/or jets, could 
have improved sensitivity. 
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